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Summary 

Absorption  spectroscopy  techniques  were  used  for  measurement  of  absolute 
concentrations  of  free  radical  species  in  premixed,  laminar,  flat-flames  of  CH4/N2O  and 
CH4/NO/O2.  Absolute  concentration  measurements  were  performed  for  the  species  OH, 
CN,  and  NH  in  63  Torr  flames.  These  absolute  concentration  measurements  were  used  to 
calibrate  relative  concentration  profiles  taken  in  previous  experiments  using  laser-induced 
fluorescence  (LIF)  spectroscopy.  Significant  absorptions  of  these  species  could  be 
obtained  even  in  the  very  short  path  length  flames  studied.  These  measurements  permit  the 
experimental  LIF  measurements  to  be  directly  compared  to  the  species  concentration 
profiles  predicted  by  the  CHEMKIN/PREMIX  flame  modeling  code. 
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Introduction 


The  measurement  of  radical  species  concentration  is  vital  in  the  study  of  combustion 
chemistry.  These  free  radicals  play  a  major  role  in  the  kinetic  mechanism  because  of  their 
high  reactivity.  Unfortunately,  due  to  this  high  reactivity,  the  very  low  concentrations  of 
these  short-lived  species  can  not  be  determined  using  conventional  intrusive  techniques 
(e.g.,  probe  sampling).  Therefore,  non-intrusive  techniques  such  as  laser-induced 
fluorescence  (LIF)  spectroscopy  have  become  valuable  tools  in  the  study  of  free  radicals  in 
combustion  and  flame  chemistry.  The  LIF  technique  allows  for  measurement  of  relative 
concentration  profiles  and  absolute  temperatures  in  the  hostile  combustion  environment. 
One  drawback  of  LIF  is  its  inability  to  readily  measure  absolute  concentrations  of  these 
species. 

A  useful  technique  for  the  measurement  of  absolute  concentrations  of  free  radicals 
in  flames  is  absorption  spectroscopy  [Ref.  1].  Absorption  spectroscopy  is  a  non-intrusive 
technique  which  does  not  perturb  the  flame  environment  By  passing  light  through  a  point 
in  the  flame  and  measuring  the  relative  intensity  and  wavelength  of  the  absorbed  light,  the 
absolute  number  of  molecules  in  the  focal  volume  can  be  determined.  By  linking  the  two 
techniques  and  matching  the  absolute  concentration  obtained  through  absorption  with  the 
corresponding  relative  values  determined  using  LIF,  the  entire  LIF  profile  can  be  converted 
into  absolute  concentrations. 

The  goal  of  this  study  was  to  develop  the  techniques  to  make  absolute  concentration 
measurements  using  absorption  spectroscopy,  and  to  perform  those  measurements  in  the 
appropriate  flame  conditions  to  support  the  current  LIF  work  being  performed  at  the 
laboratory.  The  flames  studied  in  this  experiment  were  the  same  flames  studied  in  the  LIF 
experiments:  CH4/N2O  and  CH4/NO/O2  [Ref-  2].  These  flames  were  judiciously  chosen 
because  of  their  chemical  relevancy  to  rocket  propellant  combustion  mechanisms.  By 
studying  these  “simple”  flame  systems,  the  knowledge  gained  can  be  applied  to  more 
complex  combustion  systems 

The  species  which  were  measured  in  this  experiment  were  OH,  CN,  and  NH. 
These  were  chosen  specifically  to  match  the  species  in  the  LIF  study  and  because  there  is  a 
fair  amount  of  existing  work  which  uses  these  radicals  on  which  we  could  test  and  calibrate 
our  system.  Included  in  this  report  are  the  absolute  concentrations  measured  in  the 
CH4/N2O  and  the  CH4/NO/O2  flames,  the  emission  and  absorption  profiles  measured 
during  this  experiment,  and  a  sample  calculation  for  the  determination  of  absolute 
concentrations  from  absorption  data. 
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Chamber  and  Burner 


Equipment  and  Procedures 


An  overall  schematic  of  the  experimental  apparatus  as  used  in  the  absorption  studies 
is  shown  in  Fig.  1.  The  experiment  was  conducted  at  subatmospheric  pressures  inside  the 
chamber  shown  in  Fig.  2.  The  chamber  is  evacuated  through  an  access  port  in  the  top  plate 
by  a  Welch  Duo-Seal  1397  mechanical  pump  (500  1pm  pumping  speed).  The  chamber 
pressure  is  measured  using  MKS  Baratron  capacitance  manometers  and  is  regulated  by  a 
manual  flow  value.  Optical  access  to  the  chamber  is  provided  by  four  quartz  windows. 

The  flat-flame,  premixed  burner  used  in  this  experiment  is  shown  in  Fig.  3.  The 
burner  consists  of  a  1.9  cm  stainless  steel  tube  filled  with  many  1  mm  capillary  tubes  used 
to  create  a  laminar  flow.  The  burner  is  mounted  on  a  three-axis  Aerotech  ATS  400  Series 
Positioning  Stage  Translator  assembly  controlled  by  an  Aerotech  Unidex  stage  controller. 
The  translation  stages  are  enclosed  in  a  containment  vessel  attached  to  the  bottom  of  the 
chamber  as  shown  in  Fig.  1.  The  translation  assembly  allows  1  jam  resolution  in  three 
dimensions  with  a  maximum  displacement  of  6  inches  along  each  axis.  Gases  are  supplied 
to  the  burner  through  a  stainless  steel  gas  line  which  enters  the  chamber  through  the  bottom 
of  the  containment  vessel.  The  flow  rates  of  the  reactant  gases  are  regulated  by  Vacuum 
General  DynaMass  mass  flow  controllers.  A  nitrogen  shroud  gas  is  also  admitted  to  the 
chamber  through  a  bottom  access  port  and  is  regulated  through  a  manual  value. 


Table  1.  Flame  Conditions  Employed  in  this  Study 


— aT»!<Bt4.u4iitAi.77r— 

Flame 
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II 
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1.000 

|  0.360 

atom 

The  flow  rates  of  the  gases  used  in  this  study  are  shown  in  Table  1.  The  vertical 
distance  from  the  burner  where  the  measurements  were  taken  was  chosen  to  simplify  the 
referencing  of  the  relative  profiles  to  the  absolute  concentrations.  The  relative  profiles 
obtained  using  LEF  were  examined  and  the  distance  from  the  burner  surface  to  the  location 
of  the  maximum  relative  signal  was  identified.  This  was  the  location  chosen  for  the  focal 
volume  in  the  absorption  experiment.  This  allowed  the  relative  profiles  to  be  scaled  easily 
between  zero  and  the  maximum  signal.  The  distances  from  the  burner  surface  used  for  the 
various  radicals  examined  in  this  experiment  are  shown  in  Table  2.  Also,  all  the 
experiments  were  conducted  at  63  Tore  to  match  the  conditions  of  the  LIF  studies. 


Table  2.  Distance  From  Burner  Surface 


Distance  from  Burner 

Surface  (mm) 

Radical 

Flame  I 

Flame  II  | 

OH 

CN 

NH 

■KSM 

Illumination  Source.  Optics,  and  Electronics 

The  UV  light  source  used  was  a  continuous  mode  super-quiet  xenon  lamp  with  a 
150  W  power  supply,  both  manufactured  by  Hamamatsu  Corporation.  This  combination 
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provides  continuous  UV  light  with  an  intensity  variation  of  <1%  at  a  single  wavelength 
(0.2  cm'1  bandwidth).  The  arc  lamp  was  mounted  inside  a  lamp  housing  also 
manufactured  by  Hamamatsu. 

The  point  source  xenon  arc  lamp  emission  was  collimated  using  a  precision-ground 
quartz  optical  lens.  All  lenses  used  in  this  experiment  are  precision-ground  quartz  and  have 
a  focal  length  of  25  cm.  The  lenses  were  secured  using  adjustable  optical  mounts  to  aid  in 
the  alignment  procedure.  The  mounts  were  secured  to  the  working  surface  either  with 
magnetic  mounts  or  mounting  screws.  The  continuous  wave  light  source  was  chopped 
using  a  multi-speed  chopper  from  EG&G  PARC.  The  chopped  light  was  focused  into  the 
flame  zone  and  re-collimated  after  exiting  the  chamber  using  a  pair  of  lenses.  The  image 
was  turned  horizontally  using  a  mirrored  image  rotator.  The  image  was  finally  focused 
onto  the  slit  of  the  spectrometer. 

The  spectrometer  used  in  this  experiment  is  a  Spex  1702  3/4  meter  with  a  2400 
gr/mm  grating  controlled  by  a  CD2  Compudrive.  The  spectrometer  incorporates  an  RCA 
31034A  Photomultiplier  tube  (PMT).  The  PMT  is  powered  by  an  EG&G  Ortec  556  high- 
voltage  power  supply.  The  PMT  output  along  with  the  output  from  the  chopper  were 
directed  to  a  SR510  Lock-In  Amplifier  manufactured  by  Stanford  Research  Systems.  The 
signal  output  from  the  lock-in  amplifier  was  passed  to  a  MacADIOS  A/D  card  controlled  by 
an  Apple  Macintosh  IIx  computer.  The  data  collection  rate  was  set  by  a  trigger  passed  to 
the  A/D  board  from  a  Stanford  Research  DG-535  digital  delay  generator.  Data  were 
collected  and  displayed  in  real  time  using  data  collection  software  written  in-house  and 
stored  on  disk  for  future  analysis. 
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Nitrous  Oxide  |  l^jpynamass  Flow  Controller 


When  performing  this  experiment,  care  must  be  taken  to  insure  that  the 
photomultiplier  tube  is  not  saturated.  To  insure  that  the  PMT  was  indeed  unsaturated,  a 
neutral  density  filter  was  used  to  attenuate  the  incident  light  entering  the  spectrometer.  The 
resulting  drop  in  measured  intensity  was  noted.  If  the  measured  absorption  due  to  the  filter 
was  found  to  be  less  than  the  optical  density  of  the  filter,  saturation  of  the  PMT  was 
assumed  and  the  spectrometer  slits  were  adjusted  to  eliminate  this  effect. 


The  burner  surface  was  defined  as  the  vertical  location  where  the  output  of  the  PMT 
is  equal  to  the  average  of  the  maximum  and  minimum  signal  output  obtainable  at  a  specific 
wavelength.  To  determine  this  position,  the  experimental  systems  were  first  started  but  the 
burner  was  not  ignited.  The  burner  was  positioned  using  the  translation  controller  so  that 
through  visual  inspection,  the  light  beam  was  noted  to  be  well  below  the  burner.  At  this 
time,  the  burner  was  positioned  through  visual  inspection  so  that  the  beam  intersected  the 
horizontal  center  of  the  burner  side.  This  position  represented  total  blockage  of  the  light 
beam.  This  minimum  PMT  output  was  noted  and  the  burner  was  lower  until  the  beam  was 
completely  exposed  and  the  signal  leveled  off  to  a  maximum.  The  average  of  the  minimum 
and  maximum  was  calculated  and  the  burner  was  raised  until  this  value  was  obtained.  This 
position  was  defined  to  be  the  location  where  1/2  of  the  beam  is  blocked  by  the  burner 
surface.  This  method  has  proven  to  be  a  reliable  method  for  burner  surface  location 
through  visual  inspection  of  the  beam  location  on  the  burner.  The  Unidex  translation 
control  unit  was  reset  with  this  position  representing  the  origin  (0,0,0)  and  all  further 
burner  positions  were  measured  from  this  location. 


ISwliTo  I  5m  v\h  WffiTTiitKlTiuiwil  WTCfl  wimKiIpm 


Integrated  absorptions  were  measured  by  scans  recorded  across  single  rotational 
lines.  The  absorption  profiles  for  a  single  absorption  line  of  the  spectra  were  plotted 
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relative  transmission  vs.  relative  frequency  (cm*1)  and  integrated  using  a  Macintosh 
graphics  application  (Igor  from  WaveMetrics)  to  obtain  the  area  under  the  absorption  curve 
(in  units  of  cm*1).  A  typical  absorption  profile  for  the  Qi(6)  line  of  the  OH  radical  is 
shown  in  Fig.  4. 


4  6 

Relative  Frequency  (cm'1) 


Figure  4.  Typical  Absorption  Profile 

This  area,  the  integrated  absorption,  Wj"j«,  was  related  to  the  radical  concentration  by  the 
foUowing  method  discussed  in  detail  by  Lucht  et  al.  [Ref.  1]. 

As  radiation  passes  though  the  flame,  the  reduction  in  incident  intensity  Iw  is  related 
to  the  path  length,  1,  by  Beer's  law: 

Tm  =  isaiil  =  e'k™1  (l) 

Ira(0) 

where  T^  is  the  spectral  transmission,  IM(0)  and  I(o(l)  are  the  spectral  intensities  before  the 
flame  and  after  passing  through  a  path  length  1,  respectively,  and  is  the  spectral 
absorption  coefficient.  The  spectral  absorptivity  A^  is  obtained  by  rearranging  equation  1 


(2) 


Aco  =  1  -  Tw  =  W°>  -  ML 
03  03  WO) 

The  spectral  absorption  coefficient,  k(o  is  defined  by 

ku  =  Kjj  Y(to)  0) 

where  Kji  is  the  line  strength  for  transition  from  level  j  to  level  i  and  is  independent  of  to, 

and  Y(to)  is  the  line  profile,  which  accounts  for  the  line  broadening.  Combining  equations 
1-3  gives  the  following  equation  for  the  spectral  absorptivity: 

Am  =  1  -  expf-Kji  Y(cd)  1]  <4> 

When  including  simultaneous  Doppler,  collisional,  and  natural  broadening  effects,  the  line 
profile  Y(to)  becomes  the  so-called  Voigt  profile: 


Y(co)  =  -2Hn2_  V(^a) 

VtcAwd 

where  Acod  is  the  Doppler  half  width  and  the  Voigt  function  V(^,a)  is  defined  as: 


:  a)  -  a  r  exp(-y2)dy 

’  *  ].«,  a2  +  (I  -  yf 


The  integrated  absorption,  or  equivalent  width  Wj-j-  is  defined  by  equation  7: 


Wj-j' 


■r 


Aojdco 


where  0)i*  and  o>2*  are  the  limits  of  the  spectral  line. 


Combining  equations  4  and  7  and  substituting  Kj"j'  for  Kji  where  the  j  and  i  levels  are  the 
J”  and  J’  rotational  levels,  respectively,  we  obtain: 


Wj"j’ 


■r 


{ l-exp[-Kj-j’Y((o)l]  }dco 


Using  equation  5  we  then  get: 
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* 


Using  the  identity  in  equation  10,  equation  9  can  be  put  in  terms  of  one  variable, 


Notice  that  because  the  Voigt  profile  is  symmetric  about  4=0,  integration  was  performed 
over  one-half  of  the  line  and  the  result  multiplied  by  two. 


As  it  turns  out,  this  expression  (Equation  11)  cannot  be  analytically  integrated, 
therefore  numerical  methods  must  be  employed.  A  step-wise  integration  routine  was  used 
to  evaluate  the  expression.  The  first  difficulty  in  integrating  this  equation  is  the  evaluation 
of  the  Voigt  profile.  A  technique  for  the  evaluation  of  this  integral  was  published  by 
Armstrong  [Ref.  3],  and  it  was  this  method  which  was  used.  The  Voigt  profile  was 
evaluated  for  each  stepwise  position  in  this  integration  routine.  The  Voigt  profile  is 
graphed  for  several  values  of  a  in  Fig.  5.  The  numerical  value  of  the  Voigt  function  was 
then  substituted  into  the  equation  and  the  “outer”  function  was  numerically  integrated  using 
an  extended  trapezoidal  version  of  Simpson's  rule  [Ref.  4]..  Data  sets  were  calculated  for 
varying  values  of  the  broadening  parameter,  a,  and  for  different  line  strengths  Kj"j- 


(Appendix  B).  When  the  resulting  data  was  plotted  In 


V1H2A 
A©D  J 


versus 


In 


2VIn2  ^ 

Aq>d 


for  a  given  value  of  a,  the  so-called  “curves  of  growth”  were 


generated  as  shown  in  Fig.  6.  A  listing  of  the  code  used  to  produce  this  data  is  provided  in 
Appendix  A.  The  program  was  written  using  Think  C  4.0  from  Symantec  running  on  an 
Apple  Macintosh  IIx  personal  computer. 
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Figure  6.  Curves  of  Growth  Plots  for  Various  Values  of  the  Broadening 

Parameter 
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Thus,  in  the  optically  thin  limit  the  integrated  absorption  is  directly  proportional  to  both  the 
path  length  and  the  number  density  of  the  absorber.  Figure  6  demonstrates  that  deviation 

from  the  optically  thin  limit  occurs  for  W;-:-  — — 1  >  1.  For  a  typical  doppler  width  in  a 

^  AcodJ 


flame  of  0.2  cm*1  and  a  typical  measured  value  of  the  half- width  of  the  integrated 
absorption  of  2  cm*1,  this  yields  a  peak  absorption  of  24%.  Thus,  optically  thin  limit 
conditions  hold  for  peak  absorptions  of  <20%.  As  will  be  seen  later  in  this  study,  all 
absorptions  measured  here  were  <5%  and  thus  are  operating  very  close  to  the  optically  thin 
limit  The  above  “curves  of  growth”  correction  is  very  small  in  these  studies. 


Once  the  relationship  between  actual  integrated  absorption  Wj»j*  and  line  strength 
Kj"j'  is  determined,  the  number  density  at  the  particular  line  can  be  determined  from 
equation  IS. 


Njn  =  KJ"J'  2J"+1 

fy»y  Sj"j'  Tj"j’ 


(15) 


where 


mec2 

7te2 


=  1.133  x  1012  cm*1,  fy'V'  is  the  band  oscillator  strength  for  transition  from 


lower  vibrational  level  v"  to  upper  vibrational  level  v’,  Sj"j'  is  the  rotational  line  strength 
and  Tj-J'  is  a  correction  factor  used  in  calculating  band  oscillator  strengths  to  account  for 
the  increased  inter- nuclear  separation  at  higher  rotational  quantum  numbers.  Equation  16  is 
used  to  relate  the  number  density  at  a  particular  rotational  level.  Nr,  obtained  from 
equation  15,  to  the  overall  number  density,  N. 


N  =  N,.  exp(M  [Te(n)+G(v)+F(r )]) 


(16) 


where  Te(n)  is  the  electronic  term  energy,  G(v)  is  the  vibrational  term  energy  for  level  v, 
and  F(J")  is  the  rotational  term  energy  for  level  J".  Qe.  Qv»  and  Qr  are  the  electronic, 
vibrational,  and  rotational  partition  functions  and  are  given  by  equations  17,  18,  and  19, 
respectively. 


Qe  =  X  (2S+l)<j)  exp 


-hcTe(n) 

kT 


(17) 


where  S  is  the  spin  quantum  number,  <j>  equals  2  when  A*0,  and  1  when  A=0,  and  T  is  the 
temperature  at  the  point  of  interest  in  the  flame. 


-hcG(v) 
kT  . 


(18) 


(19) 
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Bv  is  the  rotational  constant.  Substituting  the  values  from  equation  17-19  into  equation  16 
yields  the  desired  result,  the  absolute  number  of  molecules  of  the  free  radical  at  the  point  of 
interest  in  the  flame.  A  sample  calculation  for  the  OH  radical  in  the  CH4/N2O  flame  can  be 
found  in  Appendix  C. 
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Results 


Emission  Spectroscopy 

The  criteria  which  were  used  for  choosing  a  rotational  line  with  which  to  perform 
absorption  measurements  were  as  follows:  the  line  must  probe  the  ground  vibrational  state 
of  the  molecule,  the  line  must  be  isolated  and  not  blended  with  any  other  lines,  and  it  must 
be  one  of  the  stronger  lines  in  order  to  improve  the  signal-to-noise  ratio.  Emission  spectra 
for  the  three  radicals  were  taken  to  identify  candidate  spectral  lines.  The  emission  spectra 
taken  for  OH,  CN,  and  NH  are  shown  in  Figs.  7,  8,  and  9,  respectively.  The  emission 
work  was  performed  in  the  same  chamber  as  the  absorption  experiments. 

OH  Emission 

The  line  used  in  the  OH  measurements  was  the  Qi(6)  transition  of  the  0,0 
vibrational  band  of  the  A2XT- X2fl  electronic  transition  which  is  located  at  308.73  nm 
[Ref.  5].  A  portion  of  the  spectrum  obtained  in  a  58  Torr  CH4/O2  flame  with  an 
equivalence  ratio  of  1.25  over  the  wavelength  range  306  to  310  nm  is  shown  in  Fig.  7. 

CN  Emission 

The  line  used  in  the  CN  measurements  was  the  unresolved  R(15)  doublet  of  the  0,0 
vibrational  band  of  the  B2X-X2X  electronic  transition  which  is  located  at  386.34  nm  [Ref. 

6] .  At  the  resolution  employed  (=0.3  cm*1)  this  doublet  could  not  be  resolved.  The 
spectrum  for  CN  obtained  in  the  N2O  flame  (Flame  II)  between  385.5  and  388.5  nm  is 
shown  in  Fig.  8. 

NH  Emission 

The  line  used  in  the  NH  measurements  was  the  R2(8)  transition  of  the  0,0 
vibrational  band  of  the  A3ri-X3X*  electronic  transition  which  is  located  at  332.69  nm  [Ref. 

7] .  The  spectrum  for  NH  obtained  in  the  N2O  flame  (Flame  II)  between  332  and  338  nm  is 
shown  in  Fig.  9.  The  R2(8)  line  is  the  middle  line  of  a  well  resolved  triplet. 

Absorption  Spectroscopy 

Absorption  scans  were  performed  over  the  applicable  spectral  transitions  in  the 
flames  studied.  The  very  small  absorptions  found  in  these  flames  necessitated  that  the 
authors  choose  a  careful  balance  of  scan  rate  and  detector  resolution.  Slow  scan  rates  allow 
greater  signal  averaging  and  thus  a  better  signal-to-noise  level.  Scans  performed  too 
slowly  may  allow  changes  in  the  flame  to  occur.  Low  resolution  scans  permit  adjacent 
spectral  lines  to  be  readily  resolved,  but  noise  increases  significantly  using  narrow  slit 
widths.  Thus  a  balance  of  these  parameters  had  to  be  chosen  to  maximize  signal  levels, 
minimize  noise  and  permit  reasonably  short  scans.  Typical  scan  times  chosen  were  =20 
minutes  for  scans  of  6  to  18  cm*1. 

OH  Absorption 

Typical  absorption  profiles  for  the  Qi(6)  transition  in  the  two  flames  are  shown  in 
Fig.  10.  The  absorption  profiles  were  graphically  integrated  using  a  data  analysis/graphics 
application  called  Igor  (WaveMetrics).  Three  trials  were  individually  integrated  for  each 
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flame  and  the  resulting  areas  were  then  averaged.  The  average  integrated  areas  for  both 
flames  are  shown  in  Table  3.  The  standard  deviation  of  these  three  trials  was  ±0.002. 

CN  Absorption 

Typical  absorption  profiles  for  the  R(15)  doublet  transition  in  the  two  flames  are 
shown  in  Fig.  11.  Three  absorption  scans  were  performed  for  each  concentration 
measurement.  Unlike  the  trials  for  OH,  the  individual  absorption  profiles  were  first 
themselves  averaged,  and  a  Gaussian  curve  fit  was  performed.  The  average  integrated  area 
of  the  Gaussian  fit  was  then  determined.  The  average  integrated  areas  for  the  two  flames 
are  shown  in  Table  3. 

NH  Absorption 

Typical  absorption  profiles  for  the  triplet  containing  the  R2(8)  transition  in  the  N2O 
flame  is  shown  in  Fig.  12  along  with  an  isolation  on  the  R2(8)  line  itself.  The  NH 
absorption  profiles  were  treated  much  the  same  as  the  CN  profiles.  Three  individual 
absorption  scans  in  the  N2O  flame  were  first  averaged,  and  a  Gaussian  curve  fit  was 
applied.  An  average  integrated  area  was  then  determined  from  this  curve  fit.  The  average 
integrated  area  is  shown  in  Table  3.  An  absorption  signal  could  not  be  detected  in  the  NO 
flame,  therefore  only  an  upper  limit  on  the  NH  concentration  could  be  estimated. 

Table  3.  Integrated  Areas  for  OH,  CN,  and  NH  Absorption  Profiles  in 
CH4/N2O  and  CH4/NO/O2  Flames  at  63  Torr 


Average  Integrated  Areas 

Species 

CH4/N2Q  Flame 

CH4/NO/O2  Flame 

OH 

0.0842  cm'1 

0.0761  cm-1 

CN 

0.00983  cm-1 

0.00573  cm*1 

NH 

0.0104  cm*1 

- 

The  absolute  concentrations  determined  for  the  three  radicals  in  both  flames  is 
shown  in  Table  4.  No  absorption  signal  could  be  detected  for  NH  in  the  CH4/N2O  flame 
and  therefore  the  value  shown  in  the  table  is  an  estimated  upper  limit,  i.e.,  ah  estimate  of 
the  maximum  concentration  of  NH  that  could  be  in  the  flame  and  still  have  the  absorption 
profile  be  undetectable. 

Table  4.  Absolute  Concentrations  of  OH,  CN,  and  NH  in  CH4/N2O  and 

CH4/NO/O2  Flames  at  63  Torr 


Peak  Mole  Fraction 

Transition 

CH4/N2O  Flame 

CH4/NO/O2  Flame 

Ql(6),  0-0,  308.73  nm 

(1.2±0.4)  x  10-2 

(9±4)  x  10-3 

R(15)  doublet,  0-0,  386.34  nm 

(5.9±2.4)  x  lO*5 

(2.9±1.2)  x  10-5 

mEEM 

R2(8),  0-0,  332.69  nm 

(4.8±2.0)  x  10-4 

<1  x  10-4 
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Conclusions 


Absorption  spectroscopy  techniques  were  used  for  measurement  of  absolute 
concentrations  of  free  radical  species  in  premixed,  laminar,  flat-flames  of  CH4/N2O  and 
CH4/NO/O2.  Absolute  concentration  measurements  were  performed  for  the  species  OH, 
CN,  and  NH  in  63  Torr  flames.  These  absolute  concentration  measurements  were  used  to 
calibrate  relative  concentration  profiles  taken  in  previous  experiments  using  laser-induced 
fluorescence  (LIF)  spectroscopy.  Significant  absorptions  of  these  species  could  be 
obtained  even  in  the  very  short  path  length  flames  studied.  These  measurements  permitted 
the  experimental  LIF  measurements  to  be  directly  compared  to  the  species  concentration 
profiles  predicted  by  the  CHEMKIN/PREMIX  flame  modeling  code. 
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Figure  7.  Emission  Spectrum  for  the  OH  radical  from  306  to  310  nm. 


18 


Figure  9.  Emission  Spectrum  for  the  NH  radical  from  332  to  338  nm. 
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Figure  10.  Absorption  Profiles  for  the  Qi(6)  transition  of  OH. 
Top:  CH4/N2O  flame;  bottom:  CH4/NO/O2  flame. 
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Appendix  A.  Listing  of  Curves  of  Growth  Integration  Routine 

The  following  code  was  written  in  Think  C  Version  4.0  on  a  Macintosh  IIx 
personal  computer.  It  has  been  left  in  ANSI  C  format  to  allow  for  greater  portability.  The 
standard  integration  was  performed  using  a  trapezoidal  variation  on  Simpson's  rule  found 
in  Numerial  Recipes  in  C  [Ref.  4J.  The  Voigt  function  was  integrated  using  a  routine 
published  by  Armstrong  [Ref.  3]  which  was  converted  from  FORTRAN  to  C. 


The  program  generates  a  specific  column  of  data  such  as  that  found  in  Appendix  B. 
The  particular  value  of  the  broadening  parameter  is  entered  as  parameter2,  which  in  the 
sample  code  is  equal  to  2.  The  range  and  step  size  of  the  values  of  the  x-axis  of  Fig.  6, 


,  fKrr  1  2 

1  n  __  i 

k  Acod  V  *  y 


are  given  by  the  variables  MINIMUMTABLE, 


MAXIMUMTABLE,  and  TABLESTEP.  The  initial  value  is  set  in  MINIMUMTABLE,  the 
final  value  is  set  in  MAXIMUMTABLE,  and  the  incremental  step  size  between  the  two 


values  is  set  in  TABLESTEP.  The  program  outputs  a  set  of  two  column  data  to  the  screen. 
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Curves  of  Growth  Integration  Routine 
Numerical  Solver 
David  A.  Alspach 

Astronautics  Laboratory (AFSC/LSCC) 
and 

Steven  Zabarnick 

University  of  Dayton  Research  Institute 


/ 


. . ***** . . . . 

Section:  includes,  defines,  prototyping,  and  global 

initialization 

Type:  initialization  section 


Date:  30  October  1990 

Description:  includes  standard  C  library's,  sets  up  define  and 

global  variables,  and  prototypes  all  function  to  be  used 
in  the  program. 

•**«*■********************#*•*#*#**#****************■***★****■*****★********/ 


♦include  <math.h> 
♦include  <stdio.h> 


♦define 

EPS 

1.0e-5 

♦define 

JMAX 

20 

♦define 

FUNC (x) 

( ( *f unc) (x) ) 

♦define 

MINIMUMTABLE 

0.0001 

/*  First  Kjj  term  in  table  */ 

♦define 

MAXIMUMTABLE 

0.0010 

/*  Final  Kjj  term  in  table  */ 

♦define 

TABLESTEP 

0.0001 

/*  Kjj  term  increment  steps*/ 

float 

func(  float 

x)  ; 

float 

qs imp (float 

(*func) () ,  float  a,  float  b) ; 

float  - 

trapzd(fioat 

(*func)(J,  float  a,  float  b,int  n)  t 

double 

voigt (double 

x,  double  yj ; 

double 

kl (double  x. 

double  y) ; 

double 

k2 (double  x. 

double  y) ; 

double 

k3 (double  x. 

double  y) ; 

double 

w[10]; 

double 

t  [  10  ]  ; 

float 

parameterl; 

/*  Variable 

definition  for  non-dimensional  Kjj  term*/ 

float 

parameter2; 

/*  Variable 

definition  for  broadening  parameter  a  */ 

Section:  main 

Type:  program  control  loop 

Date:  30  October  1990 

Description:  sets  up  the  stepping  routines  to  generate  curve  of  growth 

tables . 

*******★*★*********************★***************************★****★******/ 


main ( ) 

( 

FILE 

Point 

SFTypeList 

SFReply 

float 

parameter2*2 .  0; 


*myFile; 

thePoint; 

types; 

reply; 

y, count ; 
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print f ("Broadening  Parameter:  %f \n", parameter2) ; 
printfC'Kjj  term\twjj  term")  ; 

for (count -MIN  IMUMTABLE;  count<=MAXIM(JMTABLE;  count -count +TABLESTEP) 

< 

parameter 1 -count ; 

y-log (qsimp (func, 0,200) ) ; 

printf  (“%f\t%f\n", log ( parameter  1)  ,  y)  ; 

> 

) 

/★*************************-********************************************'* 

Section:  qsimp 

Type:  function 

Date:  30  October  1990 

Description:  Returns  the  integral  of  the  function  func  from  a  to  b. 

The  parameters  EPS  can  be  set  to  the  desired  fractional 
accuracy  and  JMAX  so  that  2AJMAX-1  is  the  maximum  allowed 
number  of  steps.  Integration  is  performed  by  Simpson's 
rule.  Note:  from  Numerical  Recipes  in  C. 

. . . . . . . . . . . 

float  qsimp (func, a, b) 
float  a,b; 
float  (*func)  ()  ; 

( 

ir.t  j  ; 

float  s ,  st ,  ost ,  os ,  trapzd  ()  ; 
void  nrerrorO; 

ost  *  os  *  -1.0e30; 
for  ( j-1; j<-JMAX; j++)  ( 

st-t  rapzd  ( f  unc,  a,  b,  j)  ; 

S“ (4 . 0*st-ost) /3 . 0; 

if  (fabs(s-os)  <  EPS*f abs (os) )  return  s; 

os«s; 

ost-st; 

> 

nrerrorC'Too  many  steps  in  routing  QSIMP") ; 

> 

/* . . . . 

Section:  trapzd 

Type:  function 

Date:  30  October  1990 

Description:  This  routine  computes  the  n'th  stage  o'  refinement  of  an 

extended  trapezoidal  rule,  func  is  input  as  a  pointer  to 
the  function  to  be  integrated  between  limits  a  and  b,  also 
input.  When  called  with  n-1,  the  routine  returns  the 
crudest  estimate  of  the  integral.  Subsequent  call  with 
n-2,3...  (in  that  sequential  order)  will  improve  the 
accuracy  by  adding  2"n-2  additional  interior  points. 

Note:  from  Numerical  Recipes  in  C. 

********★**★★*★*★*****★**★★****■***★*************★************★*********/ 

float  trapzd(func,a,b,n) 
float  a.b; 
float  (*func)  0 ; 
int  n; 

( 

float  x, i  nm, sum, del; 
static  float  s; 
static  int  it; 
inf  j; 

if  <n— 1)  ( 

it-1;  /*  it  is  the  number  of  points  to  be  added  on  the  next  call  */ 

return  (s-0.5* (b-a) * (TUNC (a) +F0NC(b) ) ) ; 

)  else  ( 
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tnm»it ; 

del- (b-a) /tnm;  /*  This  is  the  spacing  of  the  points  to  be  added  */ 
x»a*0 . 5*del; 

for  (sum-0.0,  j=l;  j<  =  it;  j  +  +  ,  x*=>del)  sura  +  »  FUNC(x); 
it  *-  2; 

s-0 . 5* ( s+ (b-a) *sum/tnm) ;  /*  This  replaces  s  by  its  refined  value  */ 

return  s; 

) 

} 

Section:  nrerror 

Type:  function 

Date:  30  October  1990 

Description:  Numerical  recipes  standard  error  handler. 

. . 


void  nrerror <error_text) 
char  error_text ( J ; 

( 

void  exit () ; 

fprintf (stderr, "Numerical  Recipes  run-time  error ... \n”) ; 
fprintf  (stderr,  ,'%s\n,,,error_text)  ; 
fprintf (stderr,  ".. .now  exiting  to  system. .. \n") ; 
exit  (1)  ; 

) 

/***••* . . . . . . . . . . . 

Section:  func 

Type:  function 

Date:  30  October  1990 

Description:  Section  which  defines  function  to  be  integrated. 


float  func (float  x) 

( 

return  l-exp(-parameterl*voigt(x,parameter2) ) ; 

/*  Note:  The  function  contains  another  function  'voigt'  which 

itself  is  an  integral  function  dependent  on  x.  This 
function  is  integrated  using  methods  described  by 
Armstrong.  */ 

) 


Section : 

Type : 

Date : 

Description : 


voigt 

function 

30  October  1990 

Function  which  controls  the  calculation  of  the  voigt 
profiles.  This  section  defines  the  constants  w  and  t 
and  decides  which  method  of  calculation  to  use  depending 
on  which  values  of  x  and  y  are  passed. 


double  voigt (double  x,  double  y) 

( 

double  v; 

w[0]-4.62243670e-l; 
w{l]”2.86675505e-l; 
w[2]»l. 090 1720 6e-l; 
w(3] -2.4810520 9e -2; 
w[4]-3.24377334e-3; 
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w [ 5 ] »2.2S333636e-4; 
w[61»7.80255648e-6; 
w(7] -1 ,08606937e-7; 
w[8]=4. 3993409 9e-10; 
w [ 9) =2.2293936Se-13; 

t [ 0 ] =0 . 245340708; 
t [ 1 ] =0 . 737473729; 
t[2J  *>1.23407622; 
t(3]»l. 7 3853771; 
t (4 ] *2 .25497400; 
t [5] *2 .78880606; 
t  ( 6]  »3 . 34785457; 
t(7]«3. 94476404; 
t [8] *4. 60368245; 
t(9]-5. 38748089; 

if  ( (y<l .0  ii  x<4 . 0)  ||  (y<l . 8/ (x+1 .0)  ii  x>4.0)  ) 

v-kl (x, y) ; 

else  if  (y<2.5  ii  x<4.0) 
v«k2  (x,y) ; 

else 

v=*k3  (x,  y)  ; 
return  v; 

) 


Section:  kl 

Type:  function 

Date:  30  October  1990 

Description:  This  section  does  the  actual  voigt  profile  calculation 

over  the  range  ycl.0. x<4;y<1.8/ (x+1) ,x>4.  It  uses  the 
Hummer-Faddeyeva/Terent 'ev  power  series  in  y  method. 

***«»«***»•» 

double  kl (double  x,  double  y) 

( 

double  y2; 
double  ul; 

double  c ( 34 ) , coef , bnol, bno2 , bn, xl, f , funct.q, yn,dn, g; 
double  dnol.d no2; 
int  i, ii; 

y2-pow(y,2) ; 

if  ( (pow(x,2) -y2)  >  70.0) 

{ 

ul-0.0; 

) 

else 

ul-exp (-pow (x, 2) +y 2 ) *cos (2 .0*x*y) ; 

if  (x  >  5.0) 

( 

dnol»- ( . 5/pow (x, 2) + . 75/pow (x, 4) +1.875/pow(x, 6) *6. 5625/pow (x, 8) + 

2  9 . 53125/pow (x,  10) +1162 . 4218/pow (x,  12) +1055 . 7  421 /pow (x, 14) ) 
dno2» (1 . 0-dnol) / (2 .0*x) ; 

) 

else 

( 

/*  Calculating  Dawsons  function  from  Hummers  Chebyshev 
coefficients  */ 
c(0|”. 1999999999972224; 
c ( 1| 1 840000000029998 ; 
c  (  2 1  =  . 1558399999965025; 
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c [ 3 1  — - . 121664  000004  3988; 
c(4] -.08 77 08 159994 0391 ; 
c [ 5 1 —  .0585141248086907; 
c(6!=.0362157301623914; 
c[7|  —  .0208497654  3  9803  6; 
c(8J=. 0111960116346270; 
c  (9)  —.562318961 67 109e-2; 
c(10]=. 2 64 8763417226 5e -2; 
c  ( 11  ]  —  .  117  32  6707S7704e-2; 
c(121-.4899519978088e-3; 
c(13]  — . 19336308  01528e-3; 
c [ 14 1 -. 72287744 6788e-4; 
c(  151 —  .25655512497  9e-4 ; 
c(16)-. 8662073 68 4 le-5; 
c[17]  —  .2787637971 9e -5; 
c (18] -. 856687362 7e-6; 
c[19I  — .  25184  33784e-6; 
c(20]-.709360221a-7^ 
c{21]  —  .191732257e-7; 
c[22]«.49801256e-8; 
ct23J  —  .12447734e-8; 
c[24]-.2997777e-9; 
c [25] — . 696450e-10; 
c(26)=.156262e-10; 
c(27J  —  ,33897e-ll; 
c[28J=.7116e-12; 
c[29]  —  .1447e-12; 
c  [  30 ] -  . 285e-13; 
c[311  —  .55e-14; 
c [32 ] - . 10e-14 ; 
c [ 33] — .2e-15; 

/*  Chenshaws  algorithm  as  given  by  Hummer  */ 

bnol-0.0; 

bno2-0 .0; 

xl-x/5.0; 

coef-4 .0*pow(xl, 2) -2 .0; 

Cor  (i-33; i>«0;i— ) 

{ 

bn-coef *bnol-bno2+c[i] ; 

bno2-bnol; 

bnol-bn; 

> 

f-xl* (bn-bno2) ; 
dnol»l .0-2 .0*x*f ; 
dno2-f ; 
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funct«y*dnol; 
if  (y  <  1.0e-8) 

( 

return  ul-1.12837917*funct; 

) 

else 

( 

q-1 .0; 
yn-y; 

for  (i-2;i<-50;i++) 

( 

dn- (x'dnol +dno2) *(-2.0)/ (double) i; 

dno2-dnol; 

dnol-dn; 

if  (fmod(i,2) 1-0.0) 

( 

q— 3-' 

yn-yn*y2; 
g-dn  *yn; 

f unct-funct +q*g; 
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) 


if  ( tabs (g/funct)  <  l.Oe-8) 

{ 

return  ul-1 . 1283791 7*funct ; 


> 


. . . 

Section : 

k2 

Type: 

function 

Date: 

30  October  1990 

Description : 

This  section  does  the  actual  voigt  profile 

calculation 

over  the  range  1.0<y<»2.5,  x<"4.  It  uses  a  20-term 
Gauss-Hermite  integration  of  K(x,y)  as  modified  by 

integration  by  parts. 

double  k2 (double  x,  double  y) 

double 

y2; 

double 

g; 

double 

r ,  s; 

int 

i; 

y2-pow(y,2) ; 

vQ 

1 

O 

o 

for  (1- 
{ 

■0;  i<10;  i++) 

r-t(il-x; 
s-t (i] +x; 

g-g+ (4 . 0*pow (t (il ,21-2.0) * (r*atan2 (r,y) +s*atan2 (s,y) -0.5*y 

> 

> 

/*»****■******* 

(log(y2+pow(r, 2) ) +log (y2+pow(s, 2) ) ) ) 

return  0 . 318309886*g; 

«***»*t**ft**0*************i.**************tt**»****t*0«*t 

Section : 

k3 

Type: 

function 

Date: 

30  October  1990 

Description : 

This  section  does  the  actual  voigt  profile 

calculation 

over  the  range  y>-2 . 5, x<4 ;y>”l . 8/ (x+1) , x>4 . 

It  uses  a 

20-term  Gauss-Hermite  integration  of  the  original  voigt 

representation  of  K(x,y). 

double  k3 (double  x,  double  y) 

( 

double  y2; 

double  g; 

int  1; 

y2-pow (y, 2) ; 
g-0.0; 

for  (i»0;i<10;i++) 

( 

g-g+ (1 .0/ (pow ( (x-t (  i) ) , 2) +y2) +1.0/ <pow ( (x+t [ i] ) , 2) +y2) ) * w [ 1 ] 

) 

return  0 . 318309886*y+g; 

) 
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Appendix  B.  Tabulated  Curves  of  Growth  Data 

The  following  data  set  is  a  sample  of  the  “curves  of  growth”  data  produced  by  the 
program  listed  in  Appendix  A.  This  data  is  presented  to  aid  future  users  of  this  code  in  the 
testing  of  their  versions  or  modifications.  The  data  presented  represents  a  wide  range  of 
values  for  a,  a=  0  to  10.  The  “a”  values  typically  used  in  the  experiment  in  the  flames 
studied  here  were  =0.02. 
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Sample  Tabulated  Curves  of  Growth  Data 


■W 

,/Wj.j.Vln  2' 
V  Awq 

-9.21034 

a=0.00 

-9.33116 

a=0.02 

-9.33122 

a=0.l0 

-9.33147 

a=1.00 

-9.33432 

a=10.00 

-9.36345 

-8.51719 

-8.63805 

-8.63811 

-8.63835 

-8.64119 

-8.67030 

-8.11173 

-8.23262 

-8.23268 

-8.23292 

-8.23573 

-8.26484 

-7.82405 

-7.94497 

-7.94503 

-7.94527 

-7.94806 

-7.97716 

-7.60090 

-7.72186 

-7.72192 

-7.72216 

-7.72493 

-7.75401 

-7.41858 

-7.53958 

-7.53963 

-7.53986 

-7.54262 

-7.57169 

-7.26443 

-7.38546 

-7.38552 

-7.38574 

-7.38848 

-7.41754 

- 7.13090 

-7.25196 

-7.25202 

-7.25224 

-7.25497 

-7.28402 

-7.01312 

-7.13422 

-7.13427 

-7.13449 

-7.13719 

-7.16623 

-6.90776 

-7.02889 

-7.02894 

-7.02916 

-7.03184 

-7.06087 

-6.21461 

-6.33610 

-6.33614 

-6.33632 

-6.33882 

-6.36774 

-5.80914 

-5.93099 

-5.93102 

-5.93115 

-5.93347 

-5.96229 

-5.52146 

-5.64366 

-5.64368 

-5.64378 

-5.64591 

-5.67462 

-5.29832 

-5.42087 

-5.42088 

-5.42094 

-5.42288 

-5.45149 

-5.11600 

-5.23890 

-5.23890 

-5.23892 

-5.24068 

-5.26919 

-4.96185 

-5.08510 

-5.08509 

-5.08507 

-5.08665 

-5.11505 

-4.82831 

-4.95192 

-4.95190 

-4.95184 

-4.95324 

-4.98153 

-4.71053 

-4.83449 

-4.83446 

-4.83436 

-4.83557 

-4.86377 

-4.60517 

-4.72949 

-4.72944 

-4.72930 

-4.73033 

-4.75842 

-3.91202 

-4.03986 

-4.03971 

-4.03919 

-4.03838 

-4.06542 

-3.50656 

-3.63792 

-3.63765 

-3.63674 

-3.63410 

-3.66010 

-3.21888 

-3.35375 

-3.35337 

-3.35208 

-3.34761 

-3.37256 

-2.99573 

-3.13411 

-3.13363 

-3.13194 

-3.12565 

-3.14956 

-2.81341 

-2.95529 

-2.95470 

-2.95263 

-2.94451 

-2.96739 

-2.65926 

-2.80463 

-2.80393 

-2.80148 

-2.79155 

-2.81338 

-2.52573 

-2.67458 

-2.67377 

-2.67094 

-2.65920 

-2.68000 

-2.40795 

-2.56028 

-2.55936 

-2.55615 

-2.54260 

-2.56236 

-2.30259 

-2.45839 

-2.45737 

-2.45377 

-2.43842 

-2.45714 

-1.60944 

-1.79959 

-1.79751 

-1.79016 

-1.75701 

-1.76544 

-1.20397 

-1.42783 

-1.42471 

-1.41370 

-1.36317 

-1.36143 

-0.91629 

-1.17320 

-1.16908 

-1.15449 

-1.08698 

-1.07519 

-0.69315 

-0.98248 

-0.97737 

-0.95930 

-0.87522 

-0.85349 

-0.51083 

-0.83195 

-0.82588 

-0.80441 

-0.70416 

-0.67261 

-0.35668 

-0.70898 

-0.70197 

-0.67719 

-0.56116 

-0.51990 

-0.22314 

-0.60603 

-0.59810 

-0.57009 

-0.43866 

-0.38781 

-0.10536 

-0.51823 

-0.50941 

-0.47825 

-0.33179 

-0.27146 

0.00000 

-0.44228 

-0.43258 

-0.39835 

-0.23723 

-0.16754 

0.69315 

-0.01375 

0.00362 

0.06467 

0.35385 

0.51137 

1.09861 

0.17289 

0.19642 

0.27860 

0.66729 

0.90279 

1  38629 

0.27753 

0.30618 

0.40550 

0.87187 

1.17662 

1.60944 

0.34527 

0.37835 

0.49206 

1.01979 

1.38609 

1.79176 

0.39357 

0.43061 

0.55683 

1.13382 

1.55492 

1.94591 

0.43040 

0.47108 

0.60848 

1.22575 

1.69577 

2.07944 

0.45985 

0.50395 

0.65155 

1.30239 

1.81617 

2.19723 

0.4842 3 

0.53156 

0.68861 

1.36794 

1.92101 

2.30259 

0.50493 

0.55536 

0.72127 

1.42515 

2.01359 

2.99573 

0.62155 

0.69865 

0.93420 

1.78045 

2.58793 

3.40120 

0.67829 

0.77793 

1.06468 

1.98070 

2.88916 

3.68888 

0.71480 

0.83469 

1.16256 

2.12185 

3.08495 

3.91202 

0.74132 

0.87990 

1.24216 

2.23110 

3.22656 

4.09435 

0.76196 

0.91801 

1.30977 

2.32025 

3.33605 

4.24850 

0.77874 

0.95131 

1.36878 

2.39555 

3.42468 

4.38203 

0.79283 

0.98109 

1.42128 

2.46073 

3.49885 

4.49981 

0.80492 

1.00818 

1.46863 

2.51818 

3.56250 

4.60517 

0.81549 

1.03313 

1.51180 

2.56953 

3.61818 
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Appendix  C.  Sample  Calculation  for  the  OH  Radical  in  a 

CH4/N2O  Flame 

This  Appendix  demonstrates  a  sample  calculation  for  converting  integrated 
absorption  profiles  obtained  from  absorption  spectroscopy  to  absolute  concentrations  of 
free  radical  species  in  flames.  All  values,  unless  otherwise  specified,  were  taken  from 
Lucht  et  al.  [Ref.  1].  The  equation  numbers  referenced  are  from  the  main  body  of  the 
present  report 

In  this  experiment,  multiple  scans  over  single  absorption  lines  were  performed. 
The  integrated  areas  for  these  individual  trials  were  calculated,  then  an  average  area  was 
determined.  It  was  this  average  area  which  was  used  in  the  calculations.  For  OH,  the 
average  integrated  absorption  for  the  Qi(6)  transition  was  experimentally  determined  to  be: 

Wj-j*  [Ql(6)l  =  0.0842  cm-1 

It  is  generally  assumed  that  the  line  broadening  parameter,  a,  is  proportional  to  pressure 
and  inversely  proportional  to  temperature: 


Lucht  et  al.  have  fit  experimental  data  to  this  function  and  have  found  that: 

a  =  600  ^  for  the  OH  radical 

where  the  units  of  P  and  T  are  atmospheres  and  Kelvin,  respectively.  All  the  flames  in  this 
study  were  run  at  63  Torr  (0.0829  atm).  The  integrated  absorption  of  OH  in  the  CH4/N2O 
flame  was  measured  at  6.0  mm  from  the  burner  surface,  where  the  temperature  is  2600  K 
[Ref.  2).  Therefore: 

a  =  0.01913 
P  =  0.0829  atm 
T  =  2600  K 
1  =  1.9  cm 


To  find  the  number  density  Ni(6)  of  the  N"  =  6,  J"  =  6.5,  rotational  level  and  the  total  OH 
number  density  Noh- 

For  the  OH2I1(v=0)  ground  electronic  state  (from  Lucht  et  al.): 

Te  =  ° 

G(v)  =  1841.9  cm’1 

For  the  Qi(6)  transition: 


Fj(6)  =  768.1  cm-1 
Sj"j’  =  12.66 
Tj"j*  =  0.969 
foo  =  0.00096 

too  =  32381  cm-1 
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(1) 


The  Doppler  half  width  A(Od  for  OH  at  2600  K  is  defined  by: 


Acoq  =  2 


2kTln2^  ^  =  o  2877  cm-l 
moHC 2 ) 


where  moH>  the  mass  of  OH,  is  2.8  x  10*26  kg.  a  is  defined  as: 


a  =  In 


-1.412 


(2) 


which  corresponds  to  the  y-axis  of  the  “curves  of  growth”  plot.  By  graphically 
interpolating  this  value  a  for  a=0.01913,  P  corresponds  to  the  x-axis  value,  or: 


P  =  In 


<Kj"j'  1  2Vln  2  ^ 

^  Acod  Vaf  J 


-1.190 


(3) 


By  rearranging  this  equation  to  solve  for  Kyy 

Kyy  =  eP  =  0.0490  cm’2.  (4) 

2  Yin  2  1 

Substituting  this  value  for  Kjnj'  into  Equation  15  from  the  text  to  solve  for  the  number 

Ijw^2 

density  for  the  Qi(6)  line,  using - =  1.133  xlO12  cm*1. 

nt2 


Ni(6)  =  Nj-  =  ^  °2  2J.>1  =  6  5Q3  x  1013  cm-3,  (5) 

n  e2  fv-V'  Sjt  Tj"j’ 

We  can  relate  this  number  to  the  overall  number  density  of  OH  by  using  Equation  16  and 
Qe,  Qv>  Or  equal  to  4, 0.360,  and  97.52,  repectively,  from  Equations  17-19. 

N0H  =  NJ"  ^^exp(^[Te(n)  +  G(v)  +  F(J")])  =  2.815  x  1015  cm*3.  (6) 

This  is  the  final  OH  number  density.  It  can  be  converted  into  mole  fraction  using  the 
following  formula: 


Mole  fractionoH  =  yoH 


_ Nqh _ -2.815  X  1015 

Total  Number  Density  7  177  Y  ini7 
at  63  Ton- and  2600  K 


0.01202. 


(7) 
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